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' We discuss the simulation of a complex dynamical system, the so-called quantum sawtooth map 

model, on a quantum computer. We show that a quantum computer can be used to efficiently extract 
relevant physical information for this model. It is possible to simulate the dynamical localization 
' of classical chaos and extract the localization length of the system with quadratic speed up with 

^ ^ respect to any known classical computation. We can also compute with algebraic speed up the 

pL^ . diffusion coefficient and the diffusion exponent both in the regimes of Brownian and anomalous 

diffusion. Finally, we show that it is possible to extract the fidelity of the quantum motion, which 
measures the stability of the system under perturbations, with exponential speed up. 
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I. INTRODUCTION 



o 

O 

I One of the main applications of computers is the simulation of dynamical models describing the evolution of complex 
■ systems. From the viewpoint of quantum computation, quantum mechanical systems play a special role. Indeed, the 
^ i' simulation of quantum many-body problems on a classical computer is a difficult task as the size of the Hilbert space 
grows exponentially with the number of particles. For instance, if we wish to simulate a chain of n spin-i particles, the 
ip5 ' size of the Hilbert space is 2". Namely, the state of this system is determined by 2" complex numbers. As observed 
^ , by Feynman in the 1980's the growth in memory requirement is only linear on a quantum computer, which is 
itself a many-body quantum system. For example, to simulate ?i spin-i particles we only need n qubits. Therefore, 
a quantum computer operating with only a few tens of qubits could outperform a classical computer. More recently, 
a few quantum efficient algorithms have been developed for various quantum systems, ranging from some many-body 
problems 0, Q to single-particle models of quantum chaos ^ |E II] - 

Any quantum algorithm has to address the problem of efficiently extracting useful information from the quantum 
' computer wave function. The result of the simulation of a quantum system is the wave function of this system, 
encoded in the n qubits of the quantum computer. The problem is that, in order to measure all iV = 2" wave function 
coefficients by means of standard polarization measurements of the n qubits, one has to repeat the quantum simulation 
a number of times exponential in the number of qubits. This procedure would spoil any quantum algorithm, even 
in the case in which such algorithm could compute the wave function with an exponential gain with respect to any 
classical computation. Nevertheless, there are some important physical questions that can be answered in an efficient 
way, and we will discuss a few examples in this paper. 

We will discuss a quantum algorithm which efficiently simulates the quantum sawtooth map, a physical model 
with rich and complex dynamics [^. This system is characterized by very different dynamical regimes, ranging 
from integrability to chaos, and from normal to anomalous diffusion; it also exhibits the phenomenon of dynamical 
localization of classical chaotic diffusion. We will show that some important physical quantities can be extracted 
efficiently by means of a quantum computer: 
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(i) the localization length of the system, which can be extracted with a quadratic speed up with respect to any 
known classical computation 7] ; 

(ii) the diffusion coefficient and the diffusion exponent, both in the regimes of normal (Brownian) and anomalous 
diffusion; in this case we obtain an algebraic speed up] 

(iii) the fidelity of quantum motion, which characterizes the stability of the system under perturbations; for this 
quantity we achieve an exponential speed up. 

The paper is organized as follows: the properties of the sawtooth map model are discussed in Sec.^] our quantum 
algorithm simulating the quantum dynamics of this model in Sec. lIIII the quantum computation of the localized regime 
and the extraction of the localization length in Sec. IIVI the quantum simulation of the phenomena of normal and 
anomalous diffusion and the computation of the diffusion coefficient and diffusion exponent in Sec. [V] the quantum 
computation of the fidelity of quantum motion in Sec. IVII our conclusions are summarized in Sec. IVIII 



II. THE SAWTOOTH MAP 



The sawtooth map is a prototype model in the studies of classical and quantum-dynamical systems and exhibits a 
rich variety of interesting physical phenomena, from complete chaos to complete integrability, normal and anomalous 
diffusion, dynamical localization, and cantori localization. Furthermore, the sawtooth map gives a good approximation 
to the motion of a particle bouncing inside a stadium billiard (which is a well-known model of classical and quantum 
chaos). 

The sawtooth map belongs to the class of periodically driven dynamical systems, governed by the Hamiltonian 

7-2 +°° 

H{ej-T)^- + v{e) s{t-jT), (1) 

j=-oo 

where (/, 9) are conjugate action-angle variables (0 < 6* < 2tt). This Hamiltonian is the sum of two terms, H{d, I; r) — 
Ho{I) + U(9;t), where Ho{I) — P/2 is just the kinetic energy of a free rotator (a particle moving on a circle 
parametrized by the coordinate 0), while 

U{9-t)^V{9)Y,5{T-jT) (2) 

j 

represents a force acting on the particle that is switched on and off instantaneously at time intervals T. Therefore, we 
say that the dynamics described by Hamiltonian is kicked. The corresponding Hamiltonian equations of motion 
are 

dH dV{9) ^ 
dH _ 

These equations can be easily integrated and one finds that the evolution from time IT~ (prior to the ^-th kick) to 
time (/ -|- 1)T^ (prior to the (/ + l)-th kick) is described by the map 

r / = / + F{9) , 

- - (4) 

{9 = 9 + TI, 

where F{9) = —dV{9)/d9 is the force acting on the particle. 

In the following, we will consider the special case V{9) = —k{9 — tt)^ /2. This map is called the sawtooth map, since 
the force F{9) = ~dV{9)/d9 — k{9 — Tr) has a sawtooth shape, with a discontinuity at 9 = 0. By rescaling I ^ J = TI, 
the classical dynamics is seen to depend only on the parameter K — kT . Indeed, in terms of the variables (J, 9) map 
Q becomes 

( J = J + K(9 --k) . 

(5) 
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The sawtooth map exhibits sensitive dependence on initial conditions, which is the distinctive feature of classical chaos: 
any small error is amplified exponentially in time. In other words, two nearby trajectories separate exponentially, 
with a rate given by the maximum Lyapunov exponent A, defined as 



lim imfMV 

where the discrete time t — t/T measures the number of map iterations and S(t) 
6J{t) and 6d{t), we differentiate map Q, obtaining 



(6) 



^/[SJ{t)]^ + [Se{t)]^. To compute 



' SJ' 


= M 


' SJ' 




■ 1 K ' 
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1 l + K 




59 



(7) 

The iteration of map lO gives 5J{t) and 59{t) as a function of 5J{0) and 59{Q) [5J{0) and 59{Q) represent a change 
of the initial conditions] . The stability matrix M has eigenvalues ^± = ^{2 + K ± \J K'^ + 4if ) , which do not depend 
on the coordinates J and 9 and are complex conjugate for —4 < < and real for K < —4 and K > 0. Thus, the 
classical motion is stable for —4 < K < and completely chaotic for K < —4 and K > 0. For K > 0, 5{t) oc 
asymptotycally in t, and therefore the maximum Lyapunov exponent is A = ln/i_|_. Similarly, we obtain A = ln|/^_| 
for K < -4. In the stable region -4 < if < 0, A = 0. 

The sawtooth map can be studied on the cylinder [J e (— oOj+oo)], or on a torus of sinite size {—nL < J < ttL, 
where L is an integer, to assure that no discontinuities are introduced in the second equation of (O when J is taken 
modulus 2ttL). Although the sawtooth map is a deterministic system, for if > and K < —4 the motion of a 
trajectory along the momentum direction is in practice indistinguishable from a random walk. Thus, one has normal 
diffusion in the action (momentum) variable and the evolution of the distribution function /(J, <) is governed by a 
Fokker-Planck equation: 



dl 
dt 



dJ\2 dJ 



The diffusion coefficient D is defined by 



D — lim 



((AJ(<))2 



(8) 



(9) 



where A J = J — (J), and (...) denotes the average over an ensemble of trajectories. If at time t = wc take a phase 
space distribution with initial momentum Jq and random phases Q < 9 < 2tt, then the solution of the Fokker-Planck 
equation ||SJ) is given by 



1 



^/2nDt 



exp 



(J -Jo)' 
2Dt 



The width ((A J(i))^) of this Gaussian distribution grows in time, according to 

{{^J{t)f)^D{K)t. 



(10) 



(11) 



For if > 1, the diffusion coefficient is well approximated by the random phase approximation, in which we assume 
that there are no correlations between the angles (phases) 9 at different times. Hence, we have 



D{K) « ((AJi)^) = 



1 

2^ 



d9{AJif = 



1 

2^ 



(12) 



where A Ji = J — J is the change in action after a single map step. For < if < 1 diffusion is slowed, due to the 
sticking of trajectories close to broken tori (known as cantori), and we have D{K) « 3.3 if (this regime is discussed 
in 0). For —4 < if < the motion is stable, the phase space has a complex structure of elliptic islands down to 
smaller and smaller scales, and one can observe anomalous diffusion, that is, ((AJ)^) (x i", with a ^ 1 (for instance, 
a — 0.57 when if — —0.1, see Fig. 0]below). The cases if — —1, —2, —3 are integrable. 

The quantum version of the sawtooth map is obtained by means of the usual quantization rules, 9^9 and 
I I — —id/ 09 (we set h = 1). The quantum evolution in one map iteration is described by a unitary operator U, 
called the Floquet operator, acting on the wave function 



ip = U Tp = exp 



ril+l)T- 

-i I dTH{9,i;T) 

IT- 



V', 



(13) 
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where H is Hamiltonian QJ. Since the potential V{0) is switched on only at discrete times IT, it is straightforward 
to obtain 

where 1 denotes the identity operator. It is important to emphasize that, while the classical sawtooth map depends 
only on the rescaled parameter K = kT, the corresponding quantum evolution (|14|l depends on k and T separately. 
The effective Planck constant is given by hcs — T. Indeed, if we consider the operator J ~ TI (J is the quantization 
of the classical rescaled action J), we have 

[§,J]=T[§J]^iT ^ih^ff. (15) 
The classical limit hcs — > is obtained by taking k ^ oo and T ^ 0, while keeping K — kT constant. 



III. QUANTUM COMPUTING OF THE QUANTUM SAWTOOTH MAP 

In the following, we describe an exponentially efficient quantum algorithm for simulation of the map H14(l It is 
based on the forward/backward quantum Fourier transform between momentum and angle bases. Such an approach 

is convenient since the operator U, introduced in Eq. H13|) . is the product of two operators, f/^, — g*fe(^-Trl)72 g^j^^j 

Ut — e^*"^^ diagonal in the 9 and / representations, respectively. This quantum algorithm requires the following 
steps for one map iteration: 

1. We apply Uk to the wave function 'ijj{9). In order to decompose the operator Uk into one- and two-qubit gates, 
we first of all write 9 in binary notation: 



(16) 



with ai G {0, 1}. Here n is the number of qubits, so that the total number of levels in the quantum sawtooth 
map is = 2". From this expansion, we obtain 



i9-n) 





--L) 
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--) 


V 2J1 


2n J 




2n / 



(17) 



that is 



7ri)2 = 47r2 E ii 
ii j2=i 



Ij^-i ® Oj^ ® Iji+i 



(18) 



where Ij is the identity operator for the qubit j and the one-qubit operators Oj^ and Oj^ act on qubits ji and 
j2, respectively. We have 



1 i, - 1 



2J 



— 1,-, 
2n ^' 



(19) 



where {^z)j denotes the Pauli operator for the qubit j. Note that the operator Oj is diagonal in the 
computational basis {|0), |1)}. We can insert (|18|) into the unitary operator Uk, obtaining the decomposition 



„ifc(e-7ri)72 ^ 



i2 



>o 



02 



^32 + 1 



(20) 

which is the product of two-qubit gates (controlled phase-shift gates), each acting non-trivially only on the 
qubits ji and j2- In the computational basis {laj^ajj) = |00), |01), |10), |11)} each two-qubit gate can be written 
as exp(i27r^fc£'jj J2), where Dj^ j^ is a diagonal matrix: 



in 









-2^(2^ -2k) 

_^(^_^) 

(^_^)(^-^) 



(21) 



Note that decomposition H2(J|) of Uk is specific to the sawtooth map. 
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2. The change from the 6 to the / representation is obtained by means of the quantum Fourier transform, which 
requires n Hadamard gates and ^n{n — 1) controlled phase-shift gates (see, e.g., 0). 

3. In the / representation, the operator Ut has essentially the same form as the operator Uk in the 9 representation, 
and therefore it can be decomposed into n? controlled phase-shift gates, similarly to Eq. pOfl . 

4. We return to the initial 6 representation by application of the inverse quantum Fourier transform. 

Thus, overall, this quantum algorithm requires 3n'^ -f n gates per map iteration (3n^ — n controlled phase-shifts and 
2n Hadamard gates). This number is to be compared with the 0(7i2") operations required by a classical computer 
to simulate one map iteration by means of a fast Fourier transform. Thus, the quantum simulation of the quantum 
sawtooth map dynamics is exponentially faster than any known classical algorithm. Note that the resources required 
to the quantum computer to simulate the evolution of the sawtooth map are only logarithmic in the system size N . 
Of course, there remains the problem of extracting useful information from the quantum-computer wave function. 
This will be discussed in the subsequent sections. 



IV. QUANTUM COMPUTING OF DYNAMICAL LOCALIZATION 

Dynamical localization is one of the most interesting phenomena that characterize the quantum behavior of classi- 
cally chaotic systems: quantum interference effects suppress chaotic diffusion in momentum, leading to exponentially 
localized wave functions. This phenomenon was first found and studied in the quantum kicked-rotator model Toj and 
has profound analogies with Anderson localization of electronic transport in disordered materials yj-]. Dynamical 
localization has been observed experimentally in the microwave ionization of Rydberg atoms [l^ and in experiments 
with cold atoms 

Dynamical localization can be studied in the sawtooth map model. In this case, map (|14|1 is studied on the cylinder 
[/ e (— oo,+cx))], which is cut-off at a finite number A'^ of levels due to the finite quantum (or classical) computer 
memory. Similarly to other models of quantum chaos, quantum interference in the sawtooth map leads to suppression 
of classical chaotic diffusion after a break time t* . For t > t*, while the classical distribution goes on diffusing, the 
quantum distribution reaches a steady state which decays exponentially over the momentum eigenbasis: 



,2 1 

Wm = \ {m\tp)\ « - exp 



2|m — mo| 



(22) 



with mo the initial value of the momentum (the index m singles out the eigenstates of /, that is, I\m) = m|m)) [T^. 
Therefore, for t > t* only ^ ((Am)^) ~ i levels are populated. 

An estimate of t* and £ can be obtained by means of the following argument 16]. The localized wave packet has 
significant projection over about £ basis states, both in the basis of the momentum eigenstates and in the basis of the 
eigenstates of the Floquet operator U defined by Eq. This operator is unitary and therefore its eigenvalues can 

be written as exp(iAi), and the so-called quasienenergies Xi are in the interval [0, 27r[. Thus, the mean level spacing 
between "significant" quasienergy eigenstates is AE w 2t:/£. The Heisenberg principle tells us that the minimum 
time required to the dynamics to resolve this energy spacing is given by 

t* w l/AE w £. (23) 

This is the break time after which the quantum feature of the dynamics reveals. Diffusion up to time t* involves a 
number of levels given by 



^{{Am)^) ^ ^D^t* ^£, (24) 

where D„i = D/T^ is the classical diffusion coefficient, measured in number of levels. The relations H23|l and H24(l 
imply 

t*^£^ D^. (25) 

Therefore, the quantum localization length £ for the average probability distribution is approximately equal to the 
classical diffusion coefficient. For the sawtooth map, 

£^Dra~ (^73)fc' . (26) 



Note that the quantum localization can take place on a finite system only if £ is smaller than the system size N . 



6 




FIG. 1: The probability distribution over the momentum basis for the sawtooth map with n — 6 qubits, k — \/3, K — \/2, 
and initial momentum mo = 0; the time average is taken in the intervals 10 < f < 20 (full curve) and 290 < t < 300 (dashed 
curve). The straight line fit, Wm oc exp(— 2|m|/^), gives a localization length £ « 12. Note that the logarithm is base ten. 



In Fig. n (taken from we show that exponential localization can already be clearly seen with n = 6 qubits. It is 
important to stress that in a quantum computer the memory capabilities grow exponentially with the number of qubits 
(the number of levels N is equal to 2"). Therefore, already with less than 40 qubits, one could make simulations 
inaccessible to today's supercomputers. Fig. ^ shows that the exponentially localized distribution, appearing at 
t w t* , is frozen in time, apart from quantum fluctuations, which we partially smooth out by averaging over a few 
map steps. The freezing of the localized distribution can be seen from comparison of the probability distributions 
taken immediately after t* (the full curve in Fig. ^ and at a much larger time t = 300 « 25t* (the dashed curve 
in the same figure). Here the localization length is £ w 12, and classical diffusion is suppressed after a break time 
t* K, (. K, Dm, in agreement with estimates H25|I -(I26 |I [the classical diffusion coefficient is Dm ~ (7r^/3)fc^ « 9.9]. This 
quantum computation up to times of the order of £ requires a number Ng « 3n^£ ^ 10^ of one- or two-qubit quantum 
gates. 

In Fig.|21 we show a quantum computation that might be performed already with a three-qubit quantum processor. 
It is possible to compare two very different regimes, namely the localized and the ergodic regime, by varying only the 
value of the quantum parameter k, while keeping the classical parameter K constant. In both cases the wave function 
is stationary (apart from quantum fluctuations), as can be seen from the comparison of the wave function patterns at 
different times. The difference between the two cases is striking. Notice that, in this example, the localization length 
£ < 1 and one can explain the results of this simulation using perturbation theory. Indeed, we have k ^ 0.35 < 1, and 
therefore we can treat the kick t/^ as a perturbation of the free-evolution operatore Ut- The case shown in Fig. |2is 
interesting since it involves only n = 3 qubits and a few tens on quantum gates. Therefore th is q uantum computation 
seems to be accessible or close to the present capabilities of NMR-based [iSlla ion-trap quantum processors. 

We now discuss how to extract the relevant information (the value of the localization length) from a quantum 
computer simulating the sawtooth-map dynamics. The localization length can be measured by running the algorithm 
repeatedly up to time t > t*. Each run is followed by a standard projective measurement on the computational 
(momentum) basis. Since the wave function at time t can be written as 

|^(i)) =^V3(m,t)|m), (27) 

m 

with |m) momentum eigenstates, such a measurement gives outcome fh with probabihty 

Wm ^ \{fh\i,{t))\^ ^ \^{rh,t)\\ (28) 

A first series of measurements would allow us to give a rough estimate of the variance ((Am)^) of the distribution 
Wm- In turn, ((Am)^) gives a first estimate of the localization length £. After this, we can store the results of the 
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FIG. 2: The probability distribution over the momentum basis for the sawtooth map with n = 3 qubits, k — K/T = KN/2-kL, 
and initial momentum mo = 0: ergodic regime ai L — 1 (full symbols) and localized regime at L = 5 (empty symbols). Circles 
(squares) represent the wave function after t — 3 {t = 50) time steps. The dashed line represents an equally weighted wave 
function. To smooth the results, we average over ten (one hundred) different values of K £ [1.4, 1.5] for the localized (ergodic) 
case. 



measurements in histogram bins of width Sm oc £ ^ ((Am)^). Finally, the localization length is extracted from a fit 
of the exponential decay of this coarse-grained distribution over the momentum basis. Elementary statistical theory 
tells us that, in this way, the localization length can be obtained with accuracy v after the order of Xjv^ computer 
runs. It is interesting to note that it is sufficient to perform a coarse-grained measurement to generate a coarse- 
grained distribution. This means that it will be sufficient to measure the most significant qubits, and ignore those 
that would give a measurement accuracy below the coarse graining bm. Thus, the number of runs and measurements 
is independent of I. 

In Fig. 121 we report a simulation of the measurement process. In the left figure we compare the exact probabilities 
given by the wave function with the result of a complete measurement of all qubits and the result of a coarse-grained 
measurement. The histograms are built from the same number of computational runs, each followed by a projective 
measurement. The coarse-grained measurement does not resolve the thinnest structures of the exact wave function. 
However, it is still possible to extract a good estimate of the localization length £ from a fit of the exponential 
decay of the probability distribution Wm- In the right figure we compare the localization lengths, extracted from the 
complete and the coarse-grained measurements, as a function of the number Nm of projective measurements. Two 
distinct behaviors are clearly distinguishable: the localization length computed from the complete measurement of all 
qubits converges slowly to the exact value for the localization length, since a large number of projective measurement 
is required in order to resolve the exponentially decaying tails. On the contrary, the coarse-grained measurements 
approaches the exact value after a much smaller number of measurements, even though the fluctuations as a function 
of the number of measurements are quite large. 

It is possible to give a better estimate of the localization length by computing the inverse participation ratio 

L^rn rn 

The inverse participation ratio determines the number of basis states significantly populated by the wave function 
and gives an estimate of the localization length of the system. We have 1 < ^ < iV, with the limiting cases ^ = 1 
and ^ = TV corresponding to a wave function delta-peaked (Wm = ^m,mo) or uniformly spread {Wm — ^/N). In 
the localized regime, ^ « £/2. We stress that the inverse participation ratio is almost insensitive to the behavior 
of exponentially small tails of the wave function. Thus, the estimate £ « 2^ is quite accurate already with a small 
number of coarse-grained mesurement (see Fig. ISJ- 

We now come to the crucial point, of estimating the gain of quantum computation of the localization length with 
respect to classical computation. First of all, we recall that it is necessary to make about t* = 0{£) map iterations to 
obtain the localized distribution, see Eq. H25|l . This is true, both for the present quantum algorithm and for classical 
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FIG. 3: Left: Simulation of a measurement experiment for the quantum sawtooth map at n = 6, = 72, T = 2iiL/N, L = 10, 
t = 50. The thick line is the exact wave function, the thin dashed (thick full) histogram represents the result of Nm runs, 
each followed by a projective measurement of all (all except two) qubits. In both cases Nm = 5 x 10"^. Right: The estimated 
localization length £ as a function of the number Nm of projective measurements. We estimate £ by fitting the probability 
decay for the complete (circles) and coarse-grained (squares) measurements. Triangles give 2^, with the inverse partecipation 
ratio computed from the coarse-grained probability distribution. The straight line is the theoretical result £ ~ 6.8, obtained 
from Eq. Ip))l . 

computation. It is reasonable to use a basis size N — 0{£) to detect localization (say, N equal to a few times the 
localization length). In such a situation, a classical computer requires 0{£^ log^) operations to extract the localization 
length, while a quantum computer would require 0(^(log^)^) elementary gates. Indeed, both classical and quantum 
computers need to perform t k t* = 0{£) = 0{N) map iterations. Therefore, the quantum computer provides a 
quadratic speed up in computing the localization length. As we shall see in Sec. I VII the quantum computation can 
provide an exponential gain (with respect to any known classical computation) in problems that require the simulation 
of dynamics up to a time t which is independent of the number of qubits. In this case, provided that we can extract 
the relevant information in a number of measurements polynomial in the number of qubits, one should compare 
0{ti}ogN)'^) elementary gates (quantum computation) with 0{tN\ogN) elementary gates (classical computation). 

V. QUANTUM COMPUTING OF BROWNIAN AND ANOMALOUS DIFFUSION 

As we have discussed in Sec. ^ the classical sawtooth map is characterized by different diffusive behaviors in the 
chaotic and semi-integrable regimes. Quantum computers could help us to study these different regimes by simulating 
the map in the deep semiclassical region — *■ 0. Let us first show that a quantum computer would be useful 
in computing the Brownian diffusion coefRcient !?,„. For this purpose, we can repeat several times the quantum 
simulation of the sawtooth map up to a given time t, ending each run with a standard projective measurement 
in the momentum basis. This allows us to compute, up to statistical errors, ((Am)^). The diffusion coefficient is 
then obtained from Eq. H24|) as Dm ~ {{^''n)'^) /t. Therefore a computation of the diffusion coefficient up to time 
t significantly involves the order of \/Dmt momentum eigenstates (other levels are only weakly populated for times 
smaller than t and can be neglected). Thus, a basis of dimension N = 0{t^^'^) is sufficient for this computation. To 
estimate the speed up of quantum computation, one should compare 0{t{\ogN)^) = 0(iV^(log A^)^) elementary gates 
(quantum computation) with with 0(iA^log A^) — 0{N^ log(A^)) elementary gates (classical computation). This gives 
an algebraic speed up. 

We note that similar computations could be done in the regime of anomalous diffusion, in which ((AJ)^) = 
T^((Arn)^) oc t", to evaluate the exponent a, a quantity of great physical interest. Such a regime is quite complex 
in the sawtooth map: Fig. 0] shows, for the classical map, the dependence of the exponent a as a function of K. As 
can be seen from this figure, the map explores subdiffusive (a < 1) and superdiffusive {a > 1) regions, up to ballistic 
diffusion (a — 2). As required by the principle of quantum to classical correspondence, the quantum sawtooth map 
follows this behavior in the deep semiclassical regime hcff 1, up to some time scale which diverges when hcs 0. 
It is important to point out that hcs drops to zero exponentially with the number of qubits (?ioff oc 1/A^ — 1/2"), 
and therefore the deep semiclassical region can be reached with a small number of qubits. For large h^ff, one can 
also study how diffusion is modified by important quantum phenomena, like quantum tunneling, localization, and 
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K 

FIG. 4: Left: Exponent a of the anomalous diffusion ({(AJ)^) oc as a function of K for the classical sawtooth map in the 
semi-integrable regime. 

quantum resonances. 

A quantum computer could help us in obtaining the exponent a of the anomalous diffusion. In this case, 
since ((Am)^) oc i", a rough estimate of the size of the basis required for the computation up to time t is 
N = 0(t"/^). Hence, we must compare 0(t(logA^)^) = 0(A^^/"(logiV)^) elementary gates (quantum computation) 
with 0{tN logN) — 0(A^("+2)/" log(A^)) elementary gates (classical computation). The speed up is again algebraic. 

VI. QUANTUM COMPUTING OF THE FIDELITY OF QUANTUM MOTION 

The simulation of quantum dynamics up to a time t which is independent of the number of qubits is useful, for 
instance, to measure dynamical correlation functions of the form 

Cit) = {^\ A^t) BiO) 1^) = (^1 (t/t)* it(o) (7* 5(0) 1^) , (30) 

where U is the time-evolution operator for the sawtooth map. Similarly, we can efficiently compute the fidelity of 
quantum motion, which is a quantity of central interest in the study of the stability of a system under perturbations 
(see, e.g., i2^i21|, 22, 23, 24, 25, 26, j27|j and references therein). The fidelity f{t) (also called the Loschmidt echo), 
measures the accuracy with which a quantum state can be recovered by inverting, at time t, the dynamics with a 
perturbed Hamiltonian. It is defined as 

f{t)^{^\iuiru'\^). (31) 

Here the wave vector \ip) evolves forward in time with the Hamiltonian H of Eq. H13|) up to time t, and then evolves 
backward in time with a perturbed Hamiltonian 7?^ ([/<: is the corresponding time-evolution operator). For instance, 
we can perturb the parameter k in the sawtooth map as follows: k k' = k + e, with e ^ k. If the evolution 
operators U and can be simulated efficiently on a quantum computer, as is the case in most physically interesting 
situations, then the fidelity of quantum motion can be evaluated with exponential speed up with respect to known 
classical computations. The same conclusion is valid for the correlation functions H30(l . 

The fidelity can be efficiently evaluated on a quantum computer, with the only requirement of an ancilla qubit, using 
the scattering circuit drawn in Fig.jS] 28, 29] . This circuit has various important applications in quantum computing, 
including quantum state tomography and quantum spectroscopy [2^ . The circuit ends up with the measurement of 
the ancilla qubit, and we have 

(a,) = Re[TT{Wp)], (ay) = lm[Tv{Wp)], (32) 

where (ctz) and (ay) are the expectation values of the Pauli spin operators az and ay for the ancilla qubit, and W is 
a unitary operator. These two expectation values can be obtained (up to statistical errors) if one runs several times 
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FIG. 5: Scattering circuit. The top line denotes the ancilla qubit, 
and W a unitary transformation. 



the bottom hne a set of n qubits, H the Hadamard gate, 



the scattering circuit. If we set p = \il}){'4>\ and W = ?7*, it is easy to see that f{t) = |Tr(M^p)p = {(^z)'^ + (fy)^- 

For this reason, provided the quantum algorithm which implements U is efficient, as it is the case for the quantum 
sawtooth map, the fidelity can be efficiently computed by means of the circuit described in Fig. We note that 
another possible way to efficiently measure the fidelity has been proposed in . 



VII. CONCLUSIONS 



In this paper, we have discussed relevant physical examples of efficient information extraction in the quantum 
computation of a dynamical system. We have shown that a quantum computer with a small number of qubits can 
efficiently simulate the quantum localization effects, simulate both the Brownian and anomalous diffusion in the 
deep semiclassical regime, and compute the fidelity of quantum motion. We would like to stress that the simulation 
of complex dynamical systems is accessible to the first generation of quantum computers with less than 10 qubits. 
Therefore, we believe that quantum algorithms for dynamical systems deserve further studies, since they are the ideal 
software for the first quantum processors. Furthermore, we emphasize that the quantum computation of quantities 
like dynamical localization or fidelity is a demanding testing ground for quantum computers. In the first case, we 
want to simulate dynamical localization, a purely quantum phenomena which is quite fragile in the presence of noise 
[33,1^; in the latter case, fidelity is computed as a result of a sophisticated many-qubit Ramsey-type interference 
experiment. Therefore the computation of these quantities appears to be a relevant test for quantum processors 
operating in the presence of decoherence and imperfection effects. 
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